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ABSTRACT 



We describe and present initial results of a weak lensing survey of nearby 
(z < 0.1) galaxy clusters in the Sloan Digital Sky Survey (SDSS). In this first 
study, gala xy clusters a r e sele cted from the SDS S spect roscopic galaxy cluster 
catalogs of iMiller et al.l (120051) and iBerlind et al.l ( l2006l ). We report a total of 
seven individual low redshift cluster weak lensing measurements which include: 
A2048, A1767, A2244, A1066, A2199, and two clusters specifically identified with 
the C4 algorithm. Our program of weak lensing of nearby galaxy clusters in the 
SDSS will eventually reach ~ 200 clusters, making it the largest weak lensing 
survey of individual galaxy clusters to date. 



Subject headings: galaxies: clusters: general — gravitational lensing 



1. Introduction 



The number density of galaxy clusters as function of mass and redshift is a well known 
probe of cosmological parameters (IHaiman et al.ll200ll ). Several methods now exist to cali- 
brate the cluster mass (optical dynamics, galaxy infall. X-ray, Sunyaev Zel'dovich, and weak 
or strong gravitational lensing) each of which have their own advantages and drawbacks. Of 
particular interest to future optical imaging surveys is weak lensing which uses the shape 
distortion of background galaxies induced by the gravity of a foreground cluster to measure 
the cluster mass. This lensing distortion is especially useful since it is independent of the 
dynamical state of the cluster and allows the cluster halo to be probed out to very large 
radii. Much progress has been made in cluster lensing studies since its initial detection 
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( iTyson et al.lll990l ). however surprisingly a recent com pilatioii estim ates that only ~ 150 
individual clusters have been studied with weak lensing (jDahld 120071 ). Typically these clus- 
ters are imaged deeply to obtain a sufficient lensing signal but usually only cover a limited 
area around each cluster, probing the cluster halo in the inner few Mpc's. At low redshift 
in particular weak lensing has not been well studied since this requires imaging an area of 
~ 1° or more around each cluster (jjoffre et al.ll2000l ). Studying a large sample of clusters 
with weak lensing at low redshift allows for a direct comparison to other methods of mass 
calibration which have been well studied in the low redshift regime. 



T o date the largest area imaging survey is the Sloan Digital Sky Survey (SDSS) (I York et al 
20001). Previous stu dies in the SDSS have measured weak lensing by stacking clusters 
( jSheldon et al.l l2008l ) in the redshift range 0.1 < z < 0.3. In spite of the shallow imag- 



ing, weak lensing of indi vidual clusters in the SDSS is also possible (IStebbins et al.l Il996 



Gould fc Villumsenlll994l ) provided the clusters lie at low redshift. In the low redshift limit 
{zi « Zs) the lensing shear signal is 



DiDis ^ 
It oc oc Di 



where Di and Dg are the angular diameter distances of the lens and source, and Dig is the 
angular diameter distance between the lens and source. The lensing noise is 



6-yt ~ Varea ~ Di 



(2) 



so the lensing signal-to- noise ratio is equal to a constant. Theref ore by taking advant age of 
the large imaging area provided by the SDSS and "going wide" (IStebbins et al.lll996l ) weak 
lensing can be used to probe the mass of low redshift galaxy clusters. 

Using the SD SS we recently reported a weak lensing measurement of the Coma Cluster 
(iKubo et al.l 120071 ) which is the lowest redshift cluster [z = 0.0236) ever measured with 
weak lensing. Since this study we have begun a program to measure other nearby, low 
redshift clusters in the SDSS. In this letter we present the ffist results from our campaign 
using clusters selected from two publicly available spect roscopic galaxy clus ter catalogs in 
the SDSS: the C4 catalog of iMiller et aP J2OO5I ) and the ISerlind et all tood ) catalog. 



2. Data 



2.1. Imaging and Spectroscopy 



For our study we use data from the SDSS, an 8000 d eg^ imaging and spectroscopic 
survey using a dedicated 2.5m telescope (IGunn et al.l l2006l ) at Apache Point Observatory 
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in New Mexico. Imaging is obtained in a time-delay-and-integrate (or drift scan) naode i n 
five filters ugriz (IFukugita et al.lll996l ) using the SDSS imaging c amera (jGunn et al.lll998l ). 



The astrometric cahbration of the SDSS is des cribed in deta i l in Pier et ahl (2003 ) a nd tfi e 
photometric calibration pipeline is described in lTucker et al.l (120061 ) and iHogg et al.l (120011 ) . 
T argets for the SPSS sp ectroscopic survey are selected using automated algorithms described 



m 



Strauss et al.l (120021 ) and spectra are obtained using two fiber-fed double spectrographs. 



The main spectroscopic galaxy sample is complete to a magnitude of r < 17.77. 



2.2. Source Galaxies 



In our lensing analysis we use galaxies drawn from the SDSS Data Re 



ease Six ( Adelman -McCarthy et a 



20081 ). Shape measurement is performed using the PHOTO pip eline (iLupton et al.M2 001) 
which measures the shapes of objects using adaptive moments ( Bernstein fc Jarvij2002. ). 



We correct for the effec ts of the point spread f unction (PSF) using the linear PSF correction 
algorithm described in iHirata fc Seljakl (120031 ) . Source galaxies used in our lensing analysis 
are required to be detected in each of the ugriz bands, clas sified by PHOTO as g alaxies 
(type=3), and have extinction corrected model magnitudes (IStoughton et al.l 120021 ) in the 
range 18 < r < 2 1.5. Other lensing st udies in the SDSS h ave used more soph i sticat ed star- 
galaxy classifiers (jsheldon et al. 2008), however similar to Mandelbaum et al. ( 2005 ) we use 
the PHOTO classification and restrict the sample to larger galaxies with a resolution factor 
R > 0.33 (or 1.5 x PSF). Here R is given by 



R= 1- 



rrcc 

M 

vrrr 



(3) 



where M„cc and Mf^^ are the sum of the second order moments (in the CCD row and 
column directions) of the object and PSF respectively. Our catalog is further restricted in 
that we only us e shape measurements from the r band which is the filter that typically has 
the best seeing (lAdelrnan-McCarthy et al.ll2008l ). Photometric redshifts are drawn from the 
SDSS photoz2 table (jOyaizu et al.l |2008| ) and we restrict our sample to galaxies between 
0.2 < 2;phot < 0.8 with error < 0.4. 



2.3. Cluster Sample 



We have selected low redshift clusters in the SDSS from two publicly avai lable spectro- 



scopic catalogs: (1) the C4 cluster catalog of iMiller et al.l (120051 ) and (2) th e iBerlind et al. 



(120061 ) catalog. The public C4 catalog is based on the SDSS Data Release Two (lAbazajian et al. 
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20041 ) and contains 748 clusters in the redshift range 0.02 < z < 0.17 over a 2600 deg^ re- 



gion. The C4 algorithm searches for clusters in a seven dimensional parameter space which 
includes position (R.A.,decl.), redshift, and four colors {u — g,g — r,r — i,i — z) . The use of 



galax y colors in this algorithm is found to minimize cluster projection effects (IMiller et al. 



20051 ). The public C4 catalog contains three different cluster centers which include the peak 



in the C4 density field, a luminosity-weighted mean centroid, and the position of the BCG 
(brightest cluster galaxy). 



The Berlind catalog uses the friends-of-friends algorithm of iHuchra fc Gellerl (1l982h to 
search for grou ps and clusters in the SDSS NYU Value Added Galaxy Catalog (NYU-VAGC ) 
(jBlantonll2005l ) which is equivalent to the SDSS Data Release Three (lAbazajian et al.ll2005l ). 
The Berlind sample consists of three volume limited catalogs and we use the Mr20 (absolute r 
magnitude brighter than —20) catalog to search for clusters since it extends over the broadest 
redshift range 0.015 < z < 0.1. For each system the Berlind catalog measures several cluster 
parameters including an unweighted group centroid, a mean redshift, and a richness estimate. 
For this study we select five of the richest clusters from the C4 catalog and two rich clusters 
from the Berlind catalog, described further in §3.21 



3. Weak Lensing Analysis 
3.1. Mass Model 

The lensing shear due to a cluster is given by 

S(< r) - S(r) 
It = (4) 

^crit 

where E(< r) is the average projected mass density within a ra dius r, £l(r) is the proj ected 



mass density at r, and Scrit is the critical surface mass density (lMiralda-Escuddll99lf). We 
fitted the data for each cluster to a Navarro, Frenk, & White profile (NFW) (jNavarro et al. 



19961 ) which has been found in numerical simulations to provide an excellent description of 
dark matter halos, ranging from galaxy to cluster sized haloes. From this mass profile a 
virial mass (M200) can be determined from 

SOOtt o 

M200 = Pcraoo, (5) 

where the virial radius (r2oo) is defined as the radius where the density reaches a value of 200 
times pc (the critical density of the universe at the redshift of the cluster). The virial radius 
is related to the halo concentration (c2oo) by r2oo = C2oo/rs where is the scale radius, the 
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radius where the density profile changes shape (JTv 



avarro et al.lll996l). The e xpression for the 



Wright fc BrainerdI ((20001). 



tangential shear due to an NFW halo is given in 

To measure the shear in each cluster we use an unweighted shear estimator given by 

It = (6) 

where et is the tangential ellipticty and N is the number of galaxies measured in logarith- 
mically spaced annuli relative to the center of each cluster. Our use of an unweighted 
estimator is valid since our source galaxies are in the magnitude and size range where the 
shape measurement error in each galaxy is small compared to the intrinsic shape noise. We 
scaled the tangential ellipticity to a tangential shear in each bin using the shear responsivity 
TZ = 1 — (7gj.T, whe re we assume a fixed value of ctsn = 0.37 in all of our measurements 
jmrata et al.ll2004l ). 

To fit the data to an NFW model we adopt a likelihood approach. We define a binned 
likelihood given by 



where x Xnfw are 



X 



A^bins 

n 



et 



-(X-XNFW) 



(7) 



Xnfw 



^NFW 



Here et is the mean tangential ellipticity and Ug^ is the standard deviation of the mean 
tangential ellipticity. 



3.2. Results 

In Figured] we show lilkelihood plots in M200 vs. C200 for five C4 clusters (A2048, A1767, 
A2244, C4 1003, and C4 3156) and two clusters from the Berlind catalog (A1066 and A2199). 
The corresponding value of r2oo is also shown on the right panels in the figure. The maximum 
likelihood value for each cluster is indicated as the large white dot and contours represent 
the 50% (blue), 75% (green), and 87.5% (yellow) confidence regions. For all of the clusters 
the halo concentration is not well constrained, which is in general typical of weak lensing 
measurements of clusters. For both sets of clusters we find that the shear is maximized when 
centered on the cluster BCG. For the C4 clusters we use the reported BCG center in the 
catalog but for the Berlind clusters we re-center on the cluster BCG by hand. All cluster 
centers reported in Table [1] indicate the position of the cluster BCG used in the lensing 
analysis. 



- 6 - 



Mass model parameters for our sample of clusters are summarized in Table [H The weak 
lensing masses of our cluster sample varies from log(M2oo) ~ 14.2 — 14.8 with the cluster 
redshifts ranging from z = 0.0306 — 0.0990. Because of the very wide area imaging available 
in the SDSS the shear field due to each cluster can be measured out to large separation. With 
this sample we are finding that the error in the mass is slightly reduced when extending the 
outer radius out to ~ 15h~^Mpc, and this is the outer cutoff radius used in our analysis. 
Beyond this radius the mass determination does not appear to improve, and may be subject 
to the effects of large scale structure (see §3.31) . St acked measurement o f galaxy clusters have 



measured the mean shear field out to ~ 30Mpc fISheldon et al.ll2008l ). but this is the first 
time we know of that the shear field of individual clusters have been measured out to this 
large of a separation. We briefly comment on two clusters in our sample below: 



Abell 2199 is a known supercluster which also contains the cluster Abell 2197. iRines et al 



(120021 ) measured the mass of this cluster using the galaxy infall method and obtained a cen- 
tral mass of 3.2 x lO^^h~^M0 which is within la of our measurement. Our imaging data 
is not deep enough to separate the contributions of A2199 and A2197 therefore our mass 
should be interpreted as an estimate of mass of the combined system. 

Abell 2244 has a relatively small separation on the sky (~ 30') and in redshift Az ~ 
0.01 to the cluster Abell 2245. Because of the relative proximity of these two clusters we 
also cannot separate out the individual contributions and report only an estimate of the 
combined mass for the system. In our calculation of the weak lensing mass we have assumed 
the redshift of Abell 2244 {z = 0.0990). 



3.3. Projection effects 



Because weak lensing relies on the shape distortion induced on background galaxies, 
intervening clusters along th e line of sight could cause additional error in the mass estimates 
of clusters (IHoekstral 1200 ll ). In the case of the SDSS the shape distortion due to these 
background clusters is expected to be small since the SDSS imaging is relatively shallow. 
This effect could however become n on-negligible for future imaging surveys such as the Dar k 
Energy Survey (lAnnis et al.l 120051 ) and the Joint Dark Energy Mission (jAlderingj |2005| ). 
Foreground structure correlated with the cluster itsel f could also cause a d ditional er ror as 
well as cause the mass estimates to be biased upward (IMetzler et al.lll999l ) (ICenlll997l ). Our 
sample of low- redshift clusters could in principle be used to study this effect, however a 
method to correct for this effect has not been developed so we leave this to future work. 
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4. Summary 



Here we summarize the main points of this work: 



1. In addition to the Coma Cluster (IKubo et al.l 120071 ). weak lensing measurements of 
other individual nearby {z < 0.1) clusters in the SDSS are possible. 

2. We present the first results from our survey which includes seven low redshift clusters 
selected from existing SDSS spectroscopic cluster catalogs. All of these clusters have no 
previous weak lensing measurement. 

3. With the very wide imaging area provided by the SDSS we are able to probe the shear 
field of individual clusters out to ~ 15h~^Mpc, further than other previous measurements. 

The final spectroscopi c cluster catalogs usin g these two cluster finding algorithms on 
the SDSS Data Release 7 (lAbazajian et al.ll2008l ) should identify a large number of other 
low redshift clusters for which individual weak lensing measurements are possible. We are 
also studying weak lensing on non-spectroscopic cluster cat alogs including the A bell catalog 
( lAbell et al .11 19891 ) and the Giriadi optical cluster catalog (iGirardi et al.lll998l ). Results of 
these analyses and a detailed discussion of lensing systematics will be presented in future 
papers. When complete our survey should include weak lensing mass estimates for > 200 
clusters over the entire SDSS, making it the largest study of individual clusters with weak 
lensing. This will allow for a detailed comparison of weak lensing derived cluster masses 
with masses using a variety of other probes at low redshift. 
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Table 1. Cluster Data and Lensing Virial Masses 



Name 


ID* 


R.A. 


Dccl. 


z 


log(M2oo) 






(J2000) 


(J2000) 








Abell 1767 


3011 


204.0347 


59.2064 


0.0704 


14.34 


+0.28 
-0.54 


Abell 2048 


8129 


228.8088 


4.3862 


0.0949 


14.78 


+0.22 
-0.32 


Abell 2244 


3004 


255.6771 


34.0600 


0.0990 


14.46 


+0.30 
-0.56 


C4 1003 


1003 


184.4213 


3.6558 


0.0771 


14.20 


+0.36 
-1.17 


C4 3156 


3156 


258.8017 


64.3191 


0.0950 


14.34 


+0.34 
-0.80 


Abell 1066 


12289 


159.7776 


5.2098 


0.0680 


14.78 


+0.20 
-0.30 


Abell 2199 


16089 


247.1593 


39.5512 


0.0306 


14.66 


+0.22 
-0.32 



^Cluster ID given in either the Berhnd or C4 catalog. 



Note. — The first five clusters in the table are from the C4 
catalog, the last two are from the Berlind catalog. Errors bars 
on log(M2oo) are la errors, where M200 is in units of h^^M©. 



- 11 - 



Abell 1767 Abell2048 




Fig. 1. — Weak Lensing virial mass (M200) and halo concentration (c2oo) contours for each 
cluster in our sample. On the right hand side of each panel we also indicate the corresponding 
virial radius r2oo- The solid white dot in each panel is the maximum likelihood value and 
contours represent the 50% (blue), 75% (green), and 87.5% (yellow) confidence regions. 



